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Abstract. - Mesoscopic particle based fluid models, such as dissipative particle dynamics, are usually as- 
sumed to be coarse-grained representations of an underlying microscopic fluid. A fundamental question is 
whether there exists a map from microscopic particles in these systems to the corresponding coarse-grained 
particles, such that the coarse-grained system has the same bulk and transport properties as the underlying 
system. In this letter, we investigate the coarse-graining of microscopic fluids using a Voronoi type projec- 
tion that has been suggested in several studies. The simulations show that the projection fails in defining 
coarse-grained particles that have a physically meaningful connection to the microscopic fluid. In particular, 
the Voronoi projection produces identical coarse-grained equilibrium properties when applied to systems with 
different microscopic interactions and different bulk properties. 
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Introduction. - Despite the tremendous computing power 
available today at the technological high-end of super comput- 
ing and distributed computing, detailed first-principles molecu- 
lar simulations are still limited to structures and mechanism on 
small time and length scales. For instance, molecular dynamics 
(MD), a well developed framework for molecular simulations, 
is at present capable of accurate modeling of systems up to typ- 
ically millions of atoms on a time scale of about 100 ns [1]. 
This is enough to simulate small viruses [2] or molecular mo- 
tors [3], but to model larger parts of a complex biomembrane or 
an entire organelle in a biological cell, there is a gap in time and 
length scale that will not be closed in any near future by increas- 
ing the computational resources. This limits our understanding 
of organization and dynamics at mesoscopic scales-a challenge 
that applies to biology, soft and hard matter physics [4] — and 
has resulted in large efforts invested in developing simulation 
techniques that allow for exploration of the mesoscopic regime. 

In soft matter physics, mesoscopic particle based techniques 
has gained popularity as tools for simulating complex flu- 
ids. Dissipative particle dynamics, smoothed particle dynam- 
ics, and a range of descending variations are representative ex- 
amples. A basic premise of these methods is that there exists 
a valid mesoscopic particle representation of fluid that in prin- 
ciple can be derived bottom-up, with interactions that give the 
correct hydrodynamic behavior in the large system limit. The 
aim of this letter is to reexamine the validity of this premise. 
We focus on the foundation of the class of models referred to 



as dissipative particle dynamics. 

Dissipative particle dynamics (DPD) was introduced in the 
early nineties by Hoogerbrugge and Koelman [5] as a simu- 
lation method for complex fluids. The method was motivated 
partly by the need to counter the problems of broken isotropy 
and Galilean invariance in lattice gas methods (LGA), partly 
by the need of a simpler simulation scheme for complex flu- 
ids. The remedy was found in the construction of a molecu- 
lar dynamics (MD) like scheme: A set of 'fluid particles' in- 
teracting with pairwise (and central) conservative, dissipative 
and stochastic forces. This construction ensures local conserva- 
tion of linear and angular momentum, a necessity for obtaining 
correct hydrodynamic behavior in the macroscopic limit [6,7]. 
The original version of DPD was later modified by Espanol 
and Warren [8] into the form most frequently encountered in 
the literature today. In this version the dissipative and stochas- 
tic forces fulfill a fluctuation-dissipation relation such that the 
steady state solution of the simulated system (in the continuous 
time limit) is the Gibbs canonical ensemble. 

With the advantages of being a simple off-lattice and hydro- 
dynamically correct method, DPD is often mentioned as one 
of the best approaches to model and simulate complex fluids 
on different time and length scales. As DPD has become a 
popular mesoscopic simulation tool, it has been applied to an 
extensive range of systems, from suspensions [9], polymer sys- 
tems [10], phase separations [11], and membrane and vesicle 
formation [12, 13], to modeling of red blood cells [14] and the 
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life cycle of a minimal protocell [15]. The span in applica- 
tions is reflecting the flexibility of DPD, allowing it to be eas- 
ily adapted to different types of systems, on different time and 
length scales. But the same flexibility also forces an important 
question: What is the physical meaning of a DPD particle? 

Throughout the literature the DPD particles are described in 
various ways as "fluid packets" [16], ". . . fluid regions, rather 
than individual atoms and molecules" [17], and "... not a sol- 
vent molecule but a fluid element, which represent clusters of 
solvent molecules" [18]. These are just a few samples, but rep- 
resentative of how DPD is assumed to connect to an underly- 
ing microscopic system. In summary, the consensus seems to 
be that a DPD particle represents some sort of clustered repre- 
sentation of particles in the corresponding microscopic system. 
How, and more important, if this connection can be formulated 
is not known. 

Considering this background we believe it is timely and nec- 
essary to ask if the concept of mesoscopic fluid particles has a 
sound physical interpretation. In technical terms the question 
is: Does there exist a projection that applied to a microscopic 
fluid results in DPD dynamics on the mesoscopic level? The 
existence of this type of projection is crucial for particle based 
methods to have a dynamical interpretation, where the trajec- 
tories of the coarse-grained particles have a meaning beyond 
re-creating certain equilibrium thermodynamical properties of 
the system. 

This view of DPD as being a coarse-grained representation 
is often put forward in the literature, usually by referring to the 
Mori-Zwanzig projection operator framework. There are also 
several reports on the formal derivation of DPD from projec- 
tion operators, e.g. the GENERIC framework developed by 
Espafiol [19]. However, it is harder to find studies that more 
explicitly state which projection operator to use in the deriva- 
tion of DPD. The best known attempt to work out a bottom- 
up derivation of DPD with an explicit projection was done by 
Flekk0y and Coveney [20]. By assuming the DPD particles to 
be Voronoi cells, they obtain the equations of motion for the 
mesoscopic system through a formal coarse-graining of the mi- 
croscopic dynamics. While the resulting equations are DPD- 
like they do not answer if the actual dynamics of the Voroni 
cells, as given by the projection, can be approximated by the 
standard DPD equations. Another study appeared more re- 
cently, where a projection operator approach was used to de- 
rive the equations of motion for clusters of particles [21]. In 
this derivation, however, the result relies on the assumption that 
the microscopic particles cannot move between clusters, which 
seems irreconcilable with the fluid character of the microscopic 
dynamics. 

In this letter we investigate using MD simulations how the 
equilibrium and transport properties of microscopic fluids carry 
through to the mesoscopic scale. This is done with projections 
defined through Voronoi tesselation of the microscopic dynam- 
ics. In the light of the results, we discuss the implications for 
mesoscopic particle methods. 

Coarse-graining. - To obtain a coarse-grained representa- 
tion of a microscopic particle system, it is necessary to define a 




Fig. 1: Periodic Voronoi tesselation in two dimensions defined by 
coarse-grained particles (gray). Each microscopic particle (black) 
is assigned to its closest coarse-grained particle and lies within the 
Voronoi cell of its assigned particle. The Voronoi cells vary both in 
number of microscopic particles within their borders, and their shape. 

mapping from the microscopic phase space to a coarse-grained 
level, reducing the number of degrees of freedom. There are 
many ways to choose a projection. But having in mind the 
picture of DPD as 'fluid particles', a projection defined by 
Voronoi tesselation of the simulation box is arguably a repre- 
sentation that is natural, as previously pointed out in the litera- 
ture [20,22]. 

The procedure works like this: A set of coarse-grained par- 
ticles, or clusters, are placed in the same simulation box as the 
microscopic system. The number of clusters depends on the 
level of coarse-graining. Each microscopic particle is assigned 
to the coarse-grained particle closest to itself; this defines a 
Voronoi tesselation of the simulation box. In Fig. 1 this is illus- 
trated for the two dimensional case with periodic boundaries, 
with coarse-grained particles in gray and microscopic particles 
in black. Note that the shape of a Voronoi cell can vary as well 
as the number of microscopic particles within a cell. We denote 
the average cluster size, i.e. the average number of microscopic 
particles per cluster, by N. 

With each microscopic particle assigned to a coarse-grained 
particle, we define a mapping of the microscopic dynamics to 
the coarse-grained level in terms of the microscopic variables 
mass m,, position r,, and momentum p,, where i = 1 ...n de- 
notes microscopic particle ;. The coarse-grained dynamics can 
be written as 

M k = f&CrWf&fo), 

i=i ' ;=i 

P, = £&(r,)p,/£&(r,), 

1=1 !=1 

V* = P k /M k , (1) 

where Mj, Pa, and V^ are the mass, momentum, and velocity of 
cluster Ic respectively. In general, ^ (r, ) is a non-negative func- 
tion which gives the relative contribution of particle i to cluster 
k. This formulation guarantees conservation of momentum and 
mass. In the Voronoi projection §t(r,-) is 1 if microscopic par- 
ticle ; is closest to coarse-grained particle k, otherwise 0. The 
projected mass of a coarse-grained particle is just the sum of all 
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the microscopic mass within its Voronoi cell, as the projected 
momentum is the sum of the microscopic momentums within 
the same cell. 

Simulation set up. - We investigate how the equilibrium 
and transport properties of a coarse-grained fluid reflect the 
properties of the underlying system by running simulations 
along two directions: first, we examine how the coarse-grained 
system depend on the average cluster size N. This is done us- 
ing a Lennard-Jones fluid as underlying system, simulated with 
standard Velocity Verlet MD using the shifted and truncated 
pairwise Lennard-Jones potential with reduced units 



t/L-j(r)=4 
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H(r c -r), (2) 



where r is the distance between particle pairs, r c is cutoff ra- 
dius, and H(x) is the Heaviside step function which is one for 
x > and zero elsewhere. In the Lennard-Jones simulations, 
r c = 3. 

Second, we apply an identical projection, average cluster 
size N = 10, to microscopic systems with different interaction 
potentials to find out how sensitive the projection is to the de- 
tails of the underlying system. In addition to the Lennard-Jones 
fluid, we simulate fluids with quadratic and linear potentials 

U Q (r) = 0.5ar c (\-r/r c ) 2 H(r c -r) and (3) 
Ul(t) = ar c (l-r/r e )H(r e -r). (4) 

The pre-factor a defines the magnitude of the force. Two 
quadratic potential fluids are simulated: Uq { with a = 25 and 
r c = 1 and Uq 2 with a = 100 and r t = 1. The linear potential 
fluid is simulated with a = 25 and r c = 1 . 

All microscopic systems are simulated in the microcanon- 
ical ensemble (constant volume, energy and number of parti- 
cles) using periodic boundary conditions. The parameter set- 
tings are time step T = 0.005, particle density p = 0.776 and 
temperature T = 0.861 (all in reduced units), chosen so that the 
Lennard-Jones system is in the fluid regime. This corresponds 
to the situation when DPD is usually assumed to represent the 
mesoscopic behavior. All projections are done with 200 clus- 
ters, with box size and number of microscopic particles varied 
accordingly to obtain the correct average cluster size N. The 
cluster positions are initially random, and then updated using 
the velocity in Eq. (1). The simulations are run to equilibrium 
or steady state equilibrium before any measurements are ob- 
tained. All results are averaged over four measurements from 
independent simulations. Both equilibrium and transport prop- 
erties of the coarse-grained system are measured directly from 
the resulting cluster coordinates and momentums. 

Simulation results. - For a particle system interacting 
with pairwise and central forces there is a one-to-one corre- 
spondence between the radial distribution function (RDF) and 
the pairwise potential [23]. This means that the equilibrium 
properties of a particle system at given density and tempera- 
ture are uniquely defined by the RDF. In order to check that the 
projection conserves the equilibrium properties, we compare 
the coarse-grained RDFs and the bulk moduli for the differ- 
ent microscopic and coarse-grained systems. We also measure 




Fig. 2: Radial distribution functions of coarse-grained systems. The 
curves correspond to systems with a different average number of mi- 
croscopic particles per cluster (from left to right: 3, 10, 20, 30, 40). 
The Lennard-Jones fluid is the underlying system in each case. 



the mass distribution of clusters, and relate this to the standard 
DPD model. How the transport properties of the microscopic 
system translate to the projected system is examined by com- 
paring the diffusion and viscosity at each level. 

Radial distribution functions. The RDFs of a coarse- 
grained Lennard-Jones fluid, using average cluster sizes from 
N = 3 to N = 40, are presented in Fig. 2. The RDFs are char- 
acterized by having a slow, almost linear, climb from zero out- 
wards, before smoothing off to one, the slope depending on the 
size of the cluster. In terms of interactions, this kind of behav- 
ior is typically of 'soft' interactions, such as the interactions 
typically used in standard DPD, and can be attributed to an ex- 
cluded volume effect of the Voronoi cells. 

Close to zero, the RDFs flatten out. This behavior can be at- 
tributed to individual clusters having zero mass, and this region 
must therefore be treated with care. This has been described 
previously in ref. [22]. The RDF for the underlying Lennard- 
Jones fluid is given by curve (a) in Fig. 3. 

To explore the sensitivity of the projection to the underly- 
ing system, we applied the projection with average cluster size 
N = 10 to different microscopic systems. The RDFs for the mi- 
croscopic systems and the corresponding cluster systems are all 
plotted in Fig. 3. While the microscopic systems have clearly 
distinguishable RDFs, the corresponding cluster RDFs are al- 
most identical. The discrepancy in the cluster RDFs for small r 
is an artifact of the projection as discussed above. 

Bulk modulus. For the different microscopic systems, we 
estimated the bulk moduli B (or the inverse compressibility) 
using the expression 



(5) 



where p is the pressure of the system, and the partial derivative 
is evaluated at constant temperature T. To obtain a numerical 
estimate for the derivative, the microscopic systems were sim- 
ulated with densities slightly higher and lower than p = 0.776. 
For each density the pressure was measured using the virial 
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Fig. 3: Upper figure: Radial distribution functions of different micro- 
scopic systems, a) Lennard- Jones potential t/L-j; b) Quadratic poten- 
tial Uq i ; c) Quadratic potential Uq 2 ; and d) Linear potential U^. Lower 
figure: Radial distribution functions of corresponding coarse-grained 
systems, all with average cluster size N = 10. 

expansion, and the derivative was calculated as pressure differ- 
ence over density difference. The bulk moduli for the corre- 
sponding clusters were obtained in the same way. To be able 
to apply the virial expansion in this case, we derived the ef- 
fective pairwise interactions of the clusters, estimated from the 
RDFs using a inverse Monte Carlo technique [23, 24]. The re- 
sulting bulk moduli are listed in Table 1 . The microscopic sys- 
tems have different moduli, while the coarse-grained systems 
are almost identical and seemingly independent of the underly- 
ing system. 

Diffusion. The diffusion coefficients were calculated from 
the mean square displacement of each system. Projections of 
a Lennard- Jones fluid, using different average cluster sizes N, 
produce the coefficients plotted in Fig. 4 (black squares), with 
standard deviation on the scale of the marker size. The diffu- 
sion coefficient of the underlying Lennard- Jones fluid is given 
by the diamond marker. We observe that the diffusion falls off 
with increasing N. The values we get by dividing each diffu- 
sion coefficient with the empirical factor jV~ 3 are nearly con- 
stant (gray circles and dashed line), and suggest an approximate 
power-law relation between cluster size and cluster diffusion. 

The diffusion coefficients for the different microscopic sys- 
tems and corresponding cluster systems (all with N = 10) are 
listed in Table 2, together with the ratio between the coeffi- 
cients. The variation in the ratios implies that the diffusion not 



Table 1: Bulk modulus for the different microscopic systems and cor- 
responding coarse-grained systems. The numbers in parentheses are 
the uncertainty of the last digit. 



Potential 


Bm 


Be 


C/l-j 


8.1(1) 


0.25(1) 


u Ql 


5.4(2) 


0.25(1) 


Uq 2 


2.6(1) 


0.25(1) 


U l 


15.7(1) 


0.25(1) 



only scales with N, but also with the underlying dynamics. 

Viscosity. The shear viscosity was measured using a 
Poiseuille flow method [25]. The external force used to drive 
the system to a steady state flow was applied to the microscopic 
system, not affecting any steps in the coarse-graining. Mea- 
surements were done for the cluster systems for average cluster 
sizes N = 3,20,30 and 40. As a check on the influence of the 
box size on the simulations, we also measured the viscosity of 
the Lennard- Jones fluid in each case. The resulting values are 
listed in Table 3. The viscosity is consistently higher (appr. 5- 
10%) for the clusters, but does not vary significantly with size 
(at least not more than the microscopic system does). 

Mass distribution. Because the projection allows micro- 
scopic particles to move between the clusters, the sizes of the 
individual clusters fluctuated around the mean value N. Fig. 5 
contains mass distributions measured in systems with average 
cluster sizes from N = 3 to 40, all projected from the same 
Lennard-Jones fluid. The distribution becomes broader with 
increasing N. In fact, due to the fluctuating volumes in the 
Voronoi tesselation, the standard deviation scales proportion- 
ally to N (Table 3). This in contrast to the standard DPD set up, 
where all particles have constant, and usually identical, mass. 

Discussion. - In this letter we have investigated the 
feasability of a microscopic foundation of coarse-grained par- 
ticle dynamics. Specifically, we try to answer the following 
question: does a Voronoi projection produce a coarse-graining 
that can be interpreted as mesoscopic particle dynamics with 
pairwise interactions? The motivation behind the question is 
that a class of simulation techniques, related to the so called 
dissipative particle dynamics (DPD) method, is based on the 
assumption that this, or perhaps some other, projection actually 
produce well defined particle dynamics on a mesoscale. 

The results that we find can be summarized as follows: The 
self-diffusion rate of the cluster centers obey a scaling law (em- 
pirically observed); The viscosity of the projected system is 
similar to the underlying system, and does not depend on the 
number of particles; The standard deviation of average clus- 
ter mass scales approximately linearly with the cluster mass, 
and implies that the clusters do not get a more well defined 
mass as the system grows; The cluster-cluster radial distribu- 
tion function is very similar for different microscopic systems, 
also reflected in the bulk modulus (and compressibility) of the 
system. 

For all microscopic systems, the bulk modulus is much 
higher than in the corresponding coarse-grained system. This 
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Fig. 4: Diffusion of clusters plotted against average cluster size N 
(black squares). The diffusion of each cluster divided by N° is ap- 
proximately constant (gray circles and dashed line), indicating a scal- 
ing law of the diffusion. The diffusion of the underlying Lennard- 
Jones fluid is shown as a reference (diamond). 



Fig. 6: Radial distribution functions for different realizations of the 
same coarse-grained system. The coarse-grained dynamics is resulting 
from a projection of the microscopic dynamics through a weighing 
function (see text). The curves have not reached equilibrium after half 
a million time steps. 



Table 2: Diffusion coefficients for the different microscopic systems 
and corresponding coarse-grained systems. The last row shows the 
ratio between the coefficients. The numbers in parentheses are the 
uncertainty of the last digit. 

Potential D M D c D M /D C 



U L -s 
U Ql 
Uq 2 
Ul 



0.0642(3) 
0.527(4) 
0.248(1) 
1.47(3) 



0.0168(4) 
0.083(2) 
0.058(2) 
0.200(3) 



3.82 
6.35 
4.28 
7.35 



Table 3: Viscosity values flc for cluster systems with different average 
cluster size 2V, projected from the same Lennard- Jones system. The 
Lennard- Jones viscosity values )±l-j are shown as a check for box 
size sensitivity. The numbers in parentheses are the uncertainty of the 
last digit. The last column lists standard deviation in cluster size. 



3 1.94(8) 

20 1.96(4) 

30 1.98(4) 

40 2.02(4) 



1.80(5) 
1.85(3) 
1.83(3) 
1.87(1) 



1.24 
5.79 
8.50 
10.39 




20 40 60 

Average cluster size N 



80 



Fig. 5: Histograms showing the size distribution for cluster systems 
with different average cluster size N. From left to right are the his- 
tograms for systems with N = 3, 10, 20 and 40. 



difference can be partly attributed to that we do not take the 
internal pressure of the coarse-grained particles into account. 
A more serious difficulty is that the coarse-grained bulk mod- 
ulus is approximately the same for all systems, whereas at the 
microscopic level the bulk modulus varies significantly. It is 
difficult to see how the coarse-grained dynamics could simul- 
taneously have an interaction that is consistent with the micro- 
scopic dynamics and recover the correct bulk modulus in the 
macroscopic limit. This is especially worrying when we con- 
sider that the standard approach to choosing the DPD parame- 
ters is to tune the interaction strength so that the macroscopic 
bulk modulus (or equivalently, the compressibility) is correct. 

A potential source of the problems with the Voronoi projec- 
tion is that the momenta of the cluster centers changes discon- 
tinuously when a microscopic particle moves from one cluster 
to another. In order to avoid this situation, we lifted the re- 
striction that a microscopic particle only belongs to its closest 
neighboring coarse-grained particle. This was done by replac- 
ing t,k{ r i) m Eq. (1) by a weighing function which decreases ex- 
ponentially with the distance between the cluster center and the 
microscopic particle. See ref. [20] for details. A microscopic 
particle is then distributed among several coarse-grained parti- 
cles, contributing to the mass and momentum of each coarse- 
grained particle according to the weighing functions. This re- 
sults in a soft clustering rather than the sharp Voronoi partition 
used above. In Fig. 6 we show RDFs from the same coarse- 
grained system, but with different initial conditions. After half 
a million time steps, the curves have not converged. In contrast, 
the RDF of the microscopic particles converge in less than a few 
thousand time steps after an initial temperature equilibration 
period of ten thousand steps. The slow convergence is caused 
by coarse-grained particles close together having strongly cor- 
related dynamics, staying together for long periods. While not 
completely understood, it seems to be connected to the mass 
variations in clusters. Therefore, this projection is not a practi- 
cal choice for coarse-graining. 

Another type of projection that fits the fluid particle view, 
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is to have the positions of the coarse-grained particles defined 
through K-mean clustering of the underlying particles [26]. 
The main motivation behind such a projection is that it more di- 
rectly reflects the spatial density variations in the microscopic 
system, whereas the Voronoi projection does not have such a 
direct connection to the physical distribution of the underly- 
ing particles. However, as we have shown in earlier work [26], 
using this projection is problematic as it results in a coarse- 
grained dynamics which suffers from discontinuous trajecto- 
ries, and it does not conserve the momentum of the particles. 

Our main conclusion is that the Voronoi projection has se- 
vere shortcomings as a means to define effective particles at 
a coarse-grained level from the dynamics of the microscopic 
particles. However, alternative projections fail even more dras- 
tically. It is therefore an open question, does there exist a pro- 
jection from microscopic dynamic to the mesoscopic level such 
that the mesoscopic level has a particle dynamic? We view this 
as a serious challenge to the community working with meso- 
scopic particle based simulation methods. If DPD and related 
methods are to become well established the question about 
which projection does define the coarse-graining is critical. 

There are two ways to argue around the problems pointed out 
in this letter. First, it is still possible that the Voronoi projection 
does produce a useful coarse-graining, however not in terms of 
a mesoscale point particle dynamics. The meso-particles can 
instead be interpreted as nodes in a dynamic mesh represent- 
ing local mechanical variables. Smoothed dissipative particle 
dynamics [27] may be viewed as an attempt in this direction 
even though the model still uses the particle interpretation at 
the mesoscale (which is central if local momentum conserva- 
tion is to be respected). It should be remembered that with this 
interpretation it is no longer possible to measure observables 
such as compressibility and diffusion directly from the trajec- 
tories of the meso-particles or rather the nodes in the mesh. As 
a consequence most of the machinery usually applied in con- 
nection with DPD must then be abandoned. Alternatively, as 
in [20] one abandons the idea of point particle dynamics at the 
mesoscopic level and simulates the motion of the Voroni re- 
gions. 

The second remedy is to view DPD as a pure top-down mod- 
eling approach and simply ignore that there is no clear con- 
nection to the microscopic dynamics. Rather, one then treats 
each mesoscopic particle as a thermodynamic system, as in the 
GENERIC framework [19,28], in which the equation of state 
can be directly specified [29]. With this perspective it is natu- 
ral to tune the parameters in the model to fit with macroscopic 
observables. However, we must stress that in this case it is 
unclear if the DPD model has any value as a dynamic model 
outside equilibrium. On the other hand, the lack of a clear 
connection between the scales is not unique. For example the 
classical measurement problem shows that the connection be- 
tween quantum mechanics and molecular dynamics is not ex- 
actly crystal clear. However, we think that it is unlikely that 
the relation between MD and mesoscale simulation methods 
should hide similar level of complexity, especially since the two 
levels essentially deals with the same type of representation of 
the respective system. 
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